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1 Introduction 



Impurities in correlated quantum systems have attracted considerable interest recently. In par- 
ticular in one spatial dimension exactly solvable models and powerful field theoretical methods 
have provided insights into the properties of local perturbations of ideal chain systems 
Static perturbations such as scattering potentials have a profound effect on the transport prop- 
erties of quasi one-dimensional structures such as quantum wires. Impurities with internal 
degrees of freedom, e.g. a localized magnetic moment in the Kondo problem, may be screened 
due to resonances with the electrons in the Id correlated host. 

In the framework of the Quantum Inverse Scattering Method (QISM) || the construction 
of integrable models for such systems is based on inhomogeneous vertex models constructed 
from solutions to a Yang-Baxter equation. Such inhomogeneities have first been embedded 
into periodic chains by Andrei and Johannesson for the Heisenberg model and later into 
various models including the supersymmetric t—J model of interacting electrons by various 
authors |iO|-[T3|. A direct consequence of this way of construction is the lack of backscattering 
at the impurities Consideration of such impurities in a more general field theoretical 

approach has led to the conclusion that the interactions in the integrable models are fine 
tuned to a fixed point which is unstable under renormalization flow |4]]. As a consequence the 
integrable inhomogeneities lack characteristic properties of a generic potential scatterer in a Id 
system with repulsive interactions, which has been found to drive the open chain fixed point 
leading to a vanishing of the conductivity ||. 

This can be overcome by combination of these integrable inhomogeneities with a real bound- 
ary. Again the construction of such models is possible within the QISM from solutions to the 
reflection equations (RE) |15| , |16|1 imposing consistency conditions on the possible boundary 
conditions for a given bulk system. For the t-J model the simplest such (c-number) solutions 
of the RE correspond to boundary chemical potentials and boundary magnetic fields respec- 
tively [17], 18 1 . Combining these boundary matrices with solutions of the Yang-Baxter equation 



one can derive dynamical boundary impurity models: Heisenberg models |19|,p0| and Kondo 
type impurities [|l],|2l],^2[ coupled to correlated electron systems have been studied this way. 

Here we construct the most general boundary impurities that can be realized within this 
approach by combination of the known static boundary fields for the supersymmetric t-J 
model with a dynamical impurity allowing for double occupancy of the electronic orbital at its 
site. This four state impurity alone has been studied previously for periodic chains and open 



ones with reflecting ends [0, 112], E3l] . The resulting boundary terms are characterized by the 
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boundary field and in addition by a real parameter characterizing the four dimensional typical 
representation of the graded Lie algebra gl(2\l) realized on the Hilbert space of the impurity 
and by its coupling strength to the host system which is controlled by a shift in the spectral 
parameter of the corresponding vertex. Using both boundary chemical potentials and boundary 
magnetic fields this leads to two three parametric families of boundary terms. Further models 
are obtained by application of the 'projecting method' introduced recently [p^| . Finally, we 
present the Bethe Ansatz equations determining the spectra of these impurity models. 



2 Algebraic Construction 



Following JT5|,|T6| the classification of integrable boundary conditions within the QISM is based 
on representations of two algebras T± defined in terms of reflection equations (RE). For 71 this 
equation reads: 

7Z 12 (X - 12) T_(X)1l 2 \X + (i) 7L(ji) = 71(^)ft 12 (A + fi) 71(A)ft 21 (A - //) , (2.1) 
1 2 

with 71 = T ® I and 71 = I Cg) Tlf] The algebra 71 is related to 71 by an automorphism. 
Representations of 7± determine the boundary terms in the Hamiltonian at the left (right) end 
of the chain. Since these can be chosen independently it is sufficient to consider solutions of 
(2.1) to obtain a classification of the possible boundary impurities. 

The matrix TZ in ( [2.1|) solves the quantum Yang-Baxter equation (YBE): 



ft 12 (A)^ 13 (A + fi)n 2 ' 6 (fx) = n 2 '\^n r -\x + ^)n V2 {X) 



,23 



23/ 



,13/ 



(2.2) 



where the superscripts denote the spaces in the tensor product V\ (g> V 2 <S> V 3 in which TZ^ acts 
nontrivially. 

For the t-J model this 7?.-matrix is given by: 

A .... .,. i 



(^ 12 (A))S 



■6 n S J2 + 
X + i 11 12 



X + i ili2 ' 



(2.3) 



with the graded permutation operator LT^ = 5^8^(—ly a ''- ', [a] G {0, 1} denoting the grading of 
the basis states. The c-number solutions of the RE ( ]2.1| ) corresponding to this 7?.-matrix can 

\ 



be classified |p!7[ , below we shall use the diagonal ones 

( \ \ 
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/ hX+i 
hX—i 



\ 



(2.4) 
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*For the t-J model considered here these tensor products carry a grading and we have to use a graded 
version of the QISM. For details see for example p5|. 



corresponding to a boundary chemical potential p and a boundary magnetic field h (in combi- 
nation with a chemical potential), respectively [ |T8|j . 

To construct boundary impurities carrying internal degrees of freedom we combine the 
matrices from ( p.4|) with an integrable impurity which has been considered previously in a 
periodic chain [|ll],[T2|]. In the QISM this impurity is characterized by the following £-matrix: 



[) \+m + i) a+z(| + i)^' 



^ 1 — n-f —S Qi ' 

S+ 1 - nj Qi 
\ Q\ Q\ a + 2-n ) 



C2.5) 



Here n = Y^ a =i i n ° = Z~^a c ^ c ° an< ^ ^ = 2 c a&aP c P are the electron number and spin opera- 
tors on the impurity site expressed in terms of canonical fermionic creation and annihilation 
operators. The Q a are the fermionic generators of gl(2\l) in this representation which can be 
expressed in terms of projection operators (the so called 'Hubbard Operators') X afe = \a)(b\ 
with a, 6=T,|,2,0: 

Qcr = VaTlX 0a - 2a^/aX- rT2 (2.6) 
with a = ±~ corresponding to a =|, j. 

£ 34 acts on a four- dimensional quantum space and satisfies the intertwining relation: 

TZ l2 (X - n) (£ 34 (A) ® £ 34 (/i)) = (C u (fi) ® £ 34 (A)) ft 12 (A - n) . (2.7) 

Following Sklyanin |IJ operator-valued matrix solutions of the RE are obtained by "dressing" 
a c-number solution K_(X) of the RE with £ 34 , i.e. considering the product 71(A) = £ 34 (A + 
t)K_(X) (£ 34 (-A + t))" 1 in matrix space (a shift t of the spectral parameter A is consistent with 
the intertwining relation ( |2.7| )). To construct an integrable chain with this impurity placed 
on site 1 this reasoning is iterated with the /^-operators for the t-J model, i.e. C n = TZ 0n , 
n = 2, . . . , L resulting in 

71(A) = C L (X) ■ ■ ■ C 2 {X)Cf{X + t)K.{\) {C\\-X + tj) _1 {C 2 {-X)y l ■ • ■ (C L (X)y l . (2.8) 

The integrable model is now defined through the transfer matrix 

r(A) = str [K + (A)r_(A)] . (2.9) 

str (M) = J2 a (~ l)' a ^aa is the (graded) supertrace taken in matrix space. Since the purpose 
of this paper is the classification of integrable boundary terms obtained in this class we restrict 
ourselves to the simplest case of K+(\) = 1 as a representation of the algebra T + which 



A 



corresponds to a reflecting left boundary of the chain. Then the Hamiltonian is obtained by 
differentiation of the transfer matrix with respect to the spectral parameter: 



d 

H oc i— r(A)| A =o • 

This leads to the following Hamiltonian of the quantum chain 

'z-i 



(2.10) 
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UjUj+x 1 
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+ o( n i + n i+i) 



At 2 + (a + 2) 



__n_/P,h 
2 n b > 



(2.11) 



where P projects out double occupancies on the bulk sites, and 5j, rij are the electronic spin 
and number operators on site j defined as above. The boundary terms Tif, depend on the choice 
of the boundary matrix, after a unitary transformation Ti^ is given by 



2 + 2) 
pit + 
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1 + p(a + 1) } m + <J 1 + a + p ( r + — ) }n 2 



a 



+ {l+ P a}{2^ 2 -^}+p{n 2 -2}xf 

- v / ^Tlto{xfX? T + XfX? i + / i .c.} 

- v^t 2 {xT°xf-xfx; 2 + /i.c.} , 



(2.12) 



where t = y ( [l + f ] 2 + p 2 t 2 ) and £ 2 = y ( [l + p (f - l)] 2 + p 2 t 2 ^j . Note that the repre- 
sentation of <?/(2|l) entering (|2.5|) does allow for double occupancy on the first site. 

Similarly, we obtain the following boundary operator when considering the i^-matrix: 



njh 



a 



m + (a + l)n 2 - h(l + a)n^ — h 1 1 + —J n 2 | 
- l)ni T n 2 | - ni|n 2 | + /i (1 - n 2 f) X 22 
+ /i) {xf Xj 2 + /i.e.} - Vo+l{xi° X^ + /i.e.} 
-t, fv^aX|°X{ 2 + v / aTIx| X; T + 5+^ + h.c 



(2.13) 



with t c = y(i- f ) 2 + t 2 /i 2 . 

The models constructed above can be solved using the algebraic Bethe Ansatz. Starting 
from the completely filled, fully polarized state, i.e. doubly occupied impurity site 1 and all 
other sites are occupied by a spin-f electron, we find that the spectrum of TC is determined by 



the solutions of the Bethe Ansatz equations (BAE) 

M s M c 

B h e\ {L - l) {\ k ) = \[e 2 {X h - A i )e 2 (A fc + A,) J] e^(\ k - #<)e_i(A* + ^) 
IJ e -i(^ - ^)e_i(^ + Xj) = B 9 ($ t )e a {& t + t)e a {& t - t) , (2.14) 

3=1 

where e n (x) = jz^|, M c = L + 1 — iV e , M s = L — and boundary phase shifts 

1 for If f -e 2 _ 2 (^) foriT 

S fc A) = { " and 5,(0) = { ? ; . (2.15) 

-e_!_a(A) forif* [ 1 for if* 

The energy of the corresponding Bethe state is then given by the expression 

Ms , 

E = E p b > h + 2(L-2) -J2\T-T ( 2 - 16 ) 

3=1 i + 4 

with = (4a + 8)/(4t 2 + (a + 2) 2 ) and E>* = (4a + 8) / (At 2 + (a + 2) 2 ) - 

As for the periodic and the open t-J model the zero temperature ground state and the 
low-lying charged and magnetic excitations are characterized by real solutions for the A and $ 
rapidities of the BAE. 



3 Projecting method 

Recently, it has been realized that new integrable boundary Hamiltonians may be obtained 
after fine tuning of the parameters characterizing the boundary and impurity, respectively [^4 



by projection onto an invariant subspace. An important application of this procedure is a 
Kondo spin coupled to the t-J model [^l], |22], |26| . To apply the projecting method one has 



to find a decomposition of the Hilbert space of the impurity H = Hi © H 2 and fine tune the 
parameters in 71 such that one of the following conditions is satisfied (n 12 are projectors onto 
Hi, 2 ) 

niT_n 2 = o or n 2 7ini = o . (3.1) 

Starting from K v _ we find that the decomposition H = {t, 1,0} © {2} is possible for t = 
i — 1 + ~J , the resulting Hamiltonian do not lead to new models. For 7-^ 2 (a,p) the impurity 
is a scalar one and the model reduces to an open t-J model with boundary chemical potential 
at site two. The other possible model, H^i t o(a,p), turns out to be a simple reparametrisation 
of Ht,io(a = 0,p,t). 



Choosing t — i + ^ J the condition ( |3.1|) is satisfied for the decomposition H = {2, |, j 
} © {0}. Again Ti,o(a,p) is just an open t-J model with boundary chemical potential at the 
second site. On the second subspace, however, we find a two-parametric Hamiltonian: 
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n 2 



d> {xj°xf 



Xi u X 



P + /*.*} . 



(a + b + l)(3 + 2a + 6) 

The Hamiltonian 7^2,1,1 may be constructed with aid of C 3A (a = — 1) and choosing the remaining 
two parameters as t — — |(a + &)/(a + b + 2) and p = (a + 6 + 2)/(3 + 2a + b). This £-matrix 
corresponds to the one for the dual space of the fundamental three dimensional representation 
used in |27|,|28| to construct periodic t-J models with alternating impurities and usual t-J sites. 
Choosing a = 0orfe = 0a further projection is possible: H = {2} © {f, !}. H 2 corresponds 
to a boundary chemical potential at the second site. Substituting b = in Ti.^ one 



obtains TCl J°. The resulting Hamiltonian Ttf j" can be identified with a Spin-^ Kondo impurity 



introduced in 21 



Considering only one decomposition H = {2, ]} © {0, [} satisfying ( |3.1| ) is possible. 
Choosing £ = * (f — f ) we find two new boundary Hamiltonians, namely 



n 



2,T 



/> 2 



(/i + l)(/i(a + l) - 1) 
+ (ah - l)n^n 2i + h (1 - n 2T ) Xf + 



«i + (a + 1) (n 2 - /inx-f) + n 2 -p 



h 



}+H<M^{ xfx?+h - c -} (3 - 2) 



and 



h 2 



. . 1 — ah 

n 1 + (a + l)n 2 H ; n 2 t 

h 



n u n 2T 



(3.3) 



_(h + l)(h(a + l) - 1)_ 

-V^+i{xJ°xS A + /i. c .} 

The BAE for the projected Hamiltonians coincide with the ones for the original model ( 2.14 ), 
provided that the reference state used in their construction, i.e. the state |2) in the impurity 
Hilbert space is not projected out. Hence, to obtain the spectrum of TC^o one has to use different 
BAE obtained for a suitable pseudo- vacuum. Alternatively, one may consider solutions of ( 2.14Q 
after adding the complex solutions corresponding to bound states (see e.g. Refs. |7|,p9|): 



and 



h 



(3.4) 



This results in the following set of BAE for for the Hamiltonian (|3T 

M 3 M. 



-e 3+2 /h{h)ei \\ k ) = J^e 2 (Ajfc- A i )e 2 (A fc + Aj)]^[e_i(Ajfc-^)e_i(Afc + ^) 

j^k e=i 

Ms 

JJe_i(^-Aj)e_i(^ + Aj-) = e 2a -2// l (^)e 2 +2// l ((^) , (3.5) 
0=1 

where M c = L — N e and M s = L — 1 — Np The energy of the corresponding Bethe state with 
spectral parameters {A^} and {d?} is given by 



E 



^^—^^-^-Y.jkr (3-6) 



-13 4 



4 Summary 

Starting from a particular solution of the intertwining relation ( |2.7|) built from a four-dimensional 
representation of the graded Lie algebra g/(2|l) and diagonal c-number solutions Q2.4|) of the 
RE we have constructed supersymmetric t—J models with integrable boundary impurities of 
Anderson or Kondo type, i.e. with an internal degree of freedom. Due to non-zero boundary 
potentials the resulting boundary terms break the supersymmetry of the model, the most gen- 
eral ones which can be constructed this way are given in Eqs. (|2.12|) and fl2.14|) . The presence 



of the boundary allowed for fine tuning of these potentials to project the model to a remaining 
invariant subspace. In most cases this projection led to models which had been constructed 
directly before. Eqs. ( |3.2j ) and (|3T3| ) however, are novel. 

The spectra of these models have been obtained by means of the algebraic Bethe Ansatz. 
Furthermore, having the exact dependence of the ground state energy on the parameters defin- 
ing the impurity allows for the computation of local correlations which usually are not easily 



accessible from the Bethe Ansatz solution (see Refs. [12,19] for examples). 
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